Radial coherent and intelligent states of paraxial wave equation 
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Ladder operators for the radial index of the paraxial optical modes in the cylindrical coordinates are calculated. 
The operators obey the su(l,l) algebra commutation relations. Based on this Lie algebra, we found that 
coherent modes constructed as eigenstatcs of the destruction operator or resulting from the action of the 
displacement operator on the fundamental mode are different. Some properties of these two kinds of radial 
coherent modes are studied in detail. © 2012 Optical Society of America 
OCIS codes: 070.2580, 080.3645, 260.6042. 



oscillator in 1926, when he was looking for a quantum 
wave-packet with minimum momentum-position uncer- 
tainty [14]. Later on, generalized coherent states (GCS) 
were developed for different dynamical systems and sym- 
metries. In order to define the GCS, several strategies 
have been used: minimizing the uncertainty relations 
of two conjugate observables, thus creating ideal states, 
finding the group lowering operator eigenvalue, thus cre- 
ating Intelligent States (IS), or introducing a suitable 
displacement operator acting on some reference state 
(ground state, usually), thus creating proper Coherent 
States (CS) [13,15]. 

In this paper, we studied the dynamical (or hidden) sym- 
metry of the scalar paraxial wave equation, by consider- 
ing the group algebra of the LG modes. Both IS and CS 
radial modes were then constructed starting from the al- 
gebra and some features of these modes were described in 
the position representation of the state. The most com- 
mon optical beams generated by laser sources obeys the 
scalar Paraxial Wave Equation (PWE) 



Beside Linear Momentum (LM) and Spin Angular Mo- 
mentum (SAM), photons possess transverse degrees of 
freedom: Orbital Angular Momentum (OAM) and radial 
profile distribution [1]. The light OAM received a lot of 
attentions in both classical and quantum regimes during 
last two decades [2-5]. The light OAM is characterized 
by a vortex optical phase, whose topological charge is 
given by the azimuthal index m = 0, ±1, ±2, .... On the 
contrary, the Radial Index (PJ) p = 0,1,2,..., which 
is associated to the intensity distribution of the light 
in the transverse plane, has been out of interest almost 
completely Nevertheless, these two transverse degrees 
of freedom, i.e. OAM and radial intensity distribution, 
are strictly correlated, and different OAM generators 
produce specific (and different) distributions of radial 
modes. For instance, cylindrical mode converters gen- 
erate pure Laguerre-Gauss (LG) modes starting from 
the Hermite-Gauss (HG) modes, spiral phase plates and 
pitch-fork holograms generate Hypcrgcomctric-Gaussian 
(HyGG) modes, and g-plates generate other kind of ra- 
dial HyGG modes as well [6-10]. Furthermore, some type 
of OAM projectors such as the single mode optical fibers 
widely used in the light quantum regime, affect the radial 
profile independently of the OAM value. In fact, optical 
fibers are not good OAM detectors since they are sensi- 
ble only to specific radial indices [11, 12]. 
Unlike LM and OAM photon eigenstates, the radial 
modes exhibits a more complicate structure, which ren- 
der their manipulation less obvious. In this respect, to 
overcome the difficulties inherent to the radial profile of 
paraxial beans, a deeper understanding of the underly- 
ing symmetry is advisable. Such hidden symmetry, in 
fact, can be used to unveil novel and interesting fea- 
tures of the radial modes. In particular, the radial profile 
operator and its conjugate-variable can be found from 
the symmetry of the state. Among the possible radial 
states, of particular interest is the coherent state(CS). 
The CS, in fact, is the quantum state closest to the clas- 
sical one and provides a minimum uncertainty between 
two conjugate observables [13]. Erwin Schrodinger was 
the first to derive a CS for the case of the harmonic 
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_i_ are dimensionless cylindri- 



where (p, <f>, £) and Vj 
cal coordinate and the transverse Laplacian, respec- 
tively [16]. The field at any given £-plane can be 
calculated by applying a unitary propagator U to 
the pupil wave-function, \ip)( = f(Ci OWi'y where 
f7(C; C) = ex P (f (C — O^l) ■ I n the position represen- 
tation the operator U is the Fresnel propagation ker- 

ncl K(p,<f>,C,p',c/>',C) = {pA\U{t,C)\p',4>')- Hereafter, 
without loss of generality, we consider the pupil solution 
of the PWE, and keep whole calculations in £ = plane. 
LG modes, LG P;m (p, </>, 0) := (p, 0|p,m), are a solutions 
of the scalar PWE with pupil function given by 



LGp !m (p, 0,0) 
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where m is an integer number defining the OAM eigen- 
value, p > is a not negative integer number defining the 
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Fig. 1. (Color online) Simulated spectra of different RCS 
in the basis of LG modes, (a) m = 1, \a\ = 0.2 (b) m = 3, 
\a\ = 0.2 (c) m = 3, |a| = 0.5, and (d) m = 3, \a\ = 0.9. 



radial nodes in the beam transverse plane, and I/j, m '(.) 
is the generalized Laguerre polynomial. The LG modes 
are orthogonal set of PWE's solution, i.e. (p',m'|p,m) = 
^p,p'^m,m', and carry a finite power [17]. Here, we focus 
our attention to the radial node Hilbert space keeping 
the OAM number m fixed. In order to find the dynami- 
cal (or hidden) symmetry lying behind the radial nodes 
of LG modes, we need the algebra of the node space. A 
straightforward calculation shows that the radial nodes 
ladder operators for the LG modes are given by 



P+\p,m) 
P-\p,m) 



P+\V+I,m) 
V\p-l,m). 



(2) 



where V± = Qt[^ - P 2 ), Q = ±pd p + 2p 2 -{p+\m\ + 
1) are the raising and lowering operators with eigenvalues 
P+ = s/(p+ l)(p+ 1 771 1 + 1) and V- = ^p{p + \m\), re- 
spectively. It can be easily checked that the ladder oper- 
ators obey the su(l,l) Lie algebra 



[P+,P-] = -2P [Po, P±] = ±V±, 



(3) 



where To = (2p + \m\ + 1) /2. In the equations above, p 
is the second-order differential operator 
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The states |p,m) arc eigcnstates of p with eigenvalue 
p, i.e. we have p|p,m) = p|p,m). The operators act- 
ing in the £-plane can be obtained by the correspond- 
ing operators in the pupil plane C = from = 
U((, 0)ij\f o U (C, 0), where N is any operator in the node 
Hilbert space in the pupil plane. The Casimir opera- 
tor is given by C = Pi - 1/2 (p+P- +P-P+^, where 

C|p,m) = (m 2 — l)/4|p,m) [13]. By knowing the underly- 
ing symmetry, one can try to exploit further properties 
of the radial node space. We here introduce the CS of 
the Lie group as the state resulting formally from the 



displacement of the ground state [13]. Later, we will in- 
troduce also the IS as the eigenvector of the lowering 
operator [15]. Unlike in the case of the harmonic oscilla- 
tor, in the case of su(l,l) dynamics these two states are 
not coincident, 
(i) Radial coherent state: the ground state of radial 
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Fig. 2. (Color online) Simulated propagation of two dif- 
ferent RCS: (a) m = 1, and (b) m = 2 with a = 0.6. 



nodes is |0,m), so, the CS is 

H m = e(^+- a ^-)|0,m), 



(5) 



where a — tanh (£/2) e 10 is a complex number and is 
given in terms of the two real quantities (£,$)■ The co- 
herent parameter for su(l,l) is bounded by the unit 
Poincare disk, i.e. |a| < 1. We may expand the CS in 
terms of the LG modes: 



(i-M 2 )^)£. 
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where |m|! is a normalization factor. The ra- 
dial CS (RCS) depends on the OAM value as 
well. Therefore, for each values of the OAM, 
there will be a specific set of RCS. The RCS 
is not an orthogonal set of modes: m i(a'\a) m = 



((l-H 2 )(l-M 2 )/(l-aa'*) 2 ) 
ertheless, the set is complete 
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where 



(1— |ct| 2 ) 2 ^ s ^ nc w °ight 01 measure. Therefore, the 
RCS is an over-complete set solution of the radial modes. 
The RCS is a solution of the PWE and, hence, provides 
a novel type of the paraxial modes. It can be shown 
that J2^Lo \{P> m \ a )m\ 2 = 1, so, the RCS carries a finite 
power. Fig (1) shows the projection of different RCS in 
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the basis of the first hundred p number of the LG modes. 
The position representation of RCS at the pupil, 

RCS aiTO (p,0,O) = (p,<f),0\a) m , is given by 

/ — ; — ; M+i 

/2M+1 /i — a*\ 2 - 1 (1+0) 

acs aiTB p,0,o) = 4/- rir — — e "»*-W' 

y 7r |m|! \ 1 — a ) 

Due to the presence of e lm ^, photons in the RCS mode 
\a)m carry a well-defined OAM value of mh per pho- 
ton. Propagation of two different RCSs is shown in Fig. 
(2). The mode intensity vanishes as p _2 l m at the beam 
center, so that that this mode can be generated starting 
from a Gaussian beam by simply combining a phase sin- 
gularity of charge m and parabolic transmission filter of 
order \m\. When a = the RCS reduces to a subfam- 
ily of (type-I) HyGG modes with p = (the first class 
of the so-called modified LG modes) [9]. As discussed 
in Ref. [9], the HyGG modes are over-complete set of 
modes yet non-orthogonal solutions of the PWE as well. 

(ii) Radial intelligent state: the eigenvector of lower- 
ing operator is named the IS. Here we are seeking for 
a special radial intelligent state (RIS), which obeys the 
following eigenvalues problem 

V-\v) m =V\VU- (8) 
Equation (8) can be solved in the |p,m) basis yielding 




Fig. 3. (Color online) Simulated propagation of two dif- 
ferent RISs: (a) m = 1 with 77 = 1 and (b) m = 2 with 
rj = 10, respectively. The later case has several rings 
at the pupils, which during propagation collapses to a 
doughnut shape. 



\r)) m = — 7— ... 1, V] 7= , = |p,m) (9) 
Wl ( 2 M) y/pl(p+\m\)\ 

where 77 is a free complex parameter, and /«(.) 

is the modified Bessel function of order n. Like 



the RCS modes, the RIS are a not orthogo- 
nal set of modes, since we have m'i'n'Wlm = 

J w (2VW 7 ) /\J(l\ m \ (2H)7| m |(2|77'|))<J m , m -. Nev- 
ertheless, they are a complete set since we have 

^ J d 2 7 ? X m (2|77|)/ m (2|r;|) \v) m m(v\ =% (10) 

where K n {.) is the n-order modified Bessel function of 
the second kind. In order to have an easy and clear pic- 
ture of the RIS, we may calculate its representation in 
the position space, i.e. RIS(p, <t>, 0) = {p,<p,0\q) m 



RIS,,™ (p, 0,0) = J 2 e W.-P 2 j (2^p) . 

We see that the RIS is the same as the so-called Besscl- 
Gauss beam. The Bessel-Gauss beam is then the intel- 
ligent state of radial modes and is the eigenstate of the 
lowering operator in the LG base representation. Like 
the RCS, the RIS is an eigenstate of the light OAM as 
well. Nevertheless, unlike RCS, it is not shape invariant 
under free-air propagation for general parameter 77. Fig. 
(3) shows the propagation of two RIS beams with differ- 
ent 77. We see that the intensity of one of the beams is 
fraught with a dramatic change during propagation. 

In conclusion, we have shown that the dynamical sym- 
metry of the radial node index of the scalar PWE is re- 
lated to the su(l,l) Lie algebra. Based on this algebra, 
we studied two kinds of CS associated to the radial pro- 
file. Such CSs, indeed, may be used for minimizing the 
uncertainty relation between radial index number and 
its conjugate variable. 
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